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ABSTRACT The formalism of generalized Wigner transformations developped in a 
previous paper, is applied to kinetic equations of the Lindblad type for quantum harmonic 
oscillator models. It is first applied to an oscillator coupled to an equilibrium chain of 
other oscillators having nearest-neighbour interactions. The kinetic equation is derived 
without using the so called rotating-wave approximation. Then it is shown that the classical 
limit of the corresponding phase-space equation is independent of the ordering of operators 
corresponding to the inverse of the generalized Wigner transformation, provided the latter 
is involutive. Moreover, this limit equation, which conserves the probabilistic nature of the 
distribution function and obeys an H-theorem, coincides with the kinetic equation for the 
corresponding classical system, which is derived independently and is distinct from that 
usually obtained in the litterature and not sharing the above properties. Finally the same 
formalism is applied to more general model equations used in quantum optics and it is shown 
that the above results remain unaltered. 
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1. INTRODUCTION 



In a previous paper ([1], here after called paper I), we studied the properties 
of generalized Weyl and Wigner transformations of quantum operators and 
classical phase-space functions respectively and with their help we presented a 
phase-space formulation of quantum evolution equations. In the case of open 
systems interacting with a large equilibrium bath we gave an explicit phase- 
space representation of quantum kinetic equations of the Lindblad type. As it 
is well-known, such equations are satisfactory in the sense that they conserve 
the probabilistic interpretation of the density matrix and obey an H-theorem 
([2], [3], [4] section V, [5] proposition 4.3, [6] section 4.3.). They are obtained 
in the context of various approaches to kinetic theory that employ some more 
or less systematic approximation scheme to the exact dynamics of the system 
([3]- [5], [7]- [12]). In the present paper we illustrate the formalism of paper I by 
applying it to kinetic equations of this kind, for particular harmonic oscillator 
models. 

In section 2, we first apply the general kinetic equation for an open system 
weakly interacting with a bath at equilibrium, that follows in the context 
of the formal theory presented in [5], a brief account of which is given in 
Appendix 3, to the case of a classical harmonic oscillator interacting with a 
(linear) harmonic chain, assuming nearest neighbour interactions. Moreover, 
the chain is considered to be at canonical equilibrium. We show that in the 
thermodynamic limit of an infinite chain, the interaction Hamiltonian takes a 
simple separable form. The derivation uses functional methods and leads to an 
equation of the Fokker-Planck type in the phase-space of the oscillator. 

Then in section 3 the corresponding quantum system is treated in the 
context of the same formalism, leading to a kinetic equation well-known in the 
litterature but either introduced phenomenologically or derived with the aid 
of additional approximations. Its phase-space representation via an arbitrary 
generalized Wigner transformation is accordingly given. It is also shown that 
its classical limit is unique, as long as the corresponding generalized Moyal 
bracket is a deformation of the Poisson bracket, and coincides with the classical 
equation obtained in section 2, which for the sake of completeness was derived 
there independently, thus exhibiting the consistency of our formalism, at least 
in this special case. 

Finally in section 4, a more general model equation of the Lindblad type 
for a harmonic oscillator linearly coupled to a bath of other oscillators, that has 
been considered in the litterature especially in quantum optics ([13], [14] and 
references there in) is phase-space transformed with the aid of the formalism 
of paper I. This generalizes previous results of [14] and easily implies that once 
again its classical limit is unique under the condition stated above. 



2. CLASSICAL OSCILLATOR LINEARLY COUPLED TO A 
HARMONIC CHAIN 



In this section we will study a classical harmonic oscillator model and 
in the next section our results will be related to those of the corresponding 
quantum system. The model consists of a harmonic oscillator (the subsystem 
E) coupled to an infinite set of harmonic oscillators (the reservoir R), at 
canonical equilibrium. It is assumed that the interaction is (bi)linear in the 
coordinates of E and R. The Hamiltonian reads 

H = H J: +H R + X'H I (2.1) 

H^ = ^— + -kq' 2 (2.2a) 
2M 2 y y J 



H I = J2^W (2.2c) 

fe 

where A' is a coupling parameter and the remaining symbols have an obvious 
meaning. Performing the (canonical) scale transformation 

p = M~^p' q = M^q' 



we set 



p k = m ?p k q k = miq k 



/ 77? 

H = Hx + H R + \H Ij A=«/ — A' (2-3) 

tf E = ^ + ^a 2 , ^ = A (2.4a) 
Hr = 2^ p2k+ 2^ qkkk ~ iqh Hk= m 

H I = S Te k q k q e k = (2.4c) 



k 



Actually in defining Hr we may take any positive-definite symmetric matrix 
h k i and not only a codiagonal one. This is readily seen for a finite chain by 



taking an orthogonal matrix 7 mn diagonalizing h mn and making the (canonical) 
transformation 

(q'k,P'k) (Qk,Pk) 

k k 

However for an infinite system the diagonalization procedure is more involved 
since subtleties appear due to the (presummably simultaneous) existence of a 
continuous and a point spectrum of hki- Nevertheless, if it posseses a complete 
set of eigenvectors, normal coordinates can be defined in analogy with (2.6) 
below, hence (2.8) follows and consequently kinetic equation ((2.26) below) 
remains unmodified. In view of the above discussion it is reasonable to assume 
in the present case 

h k = h% = h- k > (2.5) 

where denotes the complex conjugate. 

This model has been used extensively, particularly in the case of nearest 
neighbour interactions the various calculations being performed for a finite 
chain and then taking, the thermodynamic limit (see e.g. [15] eq(l), [16] eq(18), 
[17]- [19] and indirectly [20]). It should be noted that one often considers 
as a small parameter (Brownian motion proper) and not just the coupling 
parameter A' as in this paper. 

In case of an infinite number of nonvanishing h^s it will be assumed that 
hk — ► as | ft |— >• +oo sufficiently rapidly and as a matter of fact, they will be 
defined as the coefficients of a periodic function. Normal coordinates are then 
introduced by the Fourier series 



+ QO -j pn jig 

V q k e M & Qk = ^= ddm^ lk6 (2.6a)l 

=-oo V27T J ~* 1 



k x 



+oo 



tt(0) = P* e ~ ike * ^ = 4= r d0ir(9)e ike (2.66) 

This transformation is canonical and = ~^^ e%k6 mav be considered as a 

complete orthonormal set of vectors 

J2^(0M0')=S(9-9') (2.7a) 

k,i 



(2.76) 



where 6(9) is periodic with period 2tt and Sk,i is the Kronecker delta. 
Introducing (2.6) in (2.4) we get 

Hr = \ f d9(\n(9)\ 2 +u; 2 (9)\ ( f>(9)\ 2 ) (2.8a) 

Hi = f d9u*(9)cf)(9)q (2.86) 



— 7T 

1 1 /' 7r 

u{9) = -= ^e ine & e n = ^ d9u{9)e- md (2.9) 

n= — oo 



co 2 (9) = J2 h ^ he = u 2 (-9) > (2.10) 



fc= — oo 



because of (2.5). Note that e n is real and that it is reasonable to consider that 
e n = e_ n (cf. e.g. [15] eq.(10) for the finite case), so that u(9) is real and even. 

As we have used the basis {~^i^ etk6 }i the variables (f)(9), n(9) are complex- 
valued functions of 9. Had we chosen the basis {-7==, -7= cos k9, -7= sin k9, k € 
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A/"}, the corresponding variables would be real. Note also that by (2.9), only 
the real part contributes to Hi in (2.8b). 

In view of the above discussion, we will consider in the rest of this section 
a simple version of the Hamiltonian (2.8) namely 

H R = \j du)(TX 2 (u))+u) 2 (j) 2 (u))) (2.11a) 

Hi = J du}u(u))<j)(u})q = Wq (2.116) 

Where it(uj), (f)(iv) are real, the intergration interval is symmetric about the 
origin and u(lv) is real and even (this is not an essential assumption; see end 
of Appendix 1). It will become evident that results for the general case (2.8) 
follow by replacing u(uj) with u 2 (u))j^ in the final expressions (see also the 
discussion in Appendix 1). 

Clearly Hamilton's equations for the total system are 

q(t)=p(t) (2.12a) 
pit) = -ft 2 q -X j dwu(u))(j)(u)) (2.126) 

^ = faR = (2 ' 12c) 



tt(u;) = - — — = -uo 2 (j)(uo)-\u(uo)q (2.12d) 

using functional derivatives. Because the Hamiltonian is quadratic and the 
reservoir in a canonical equilibrium state) functional methods (see e.g. [24], [25]) 
will be used in the following to treat ab initio the infinite system. 

Because of its linearity, the model is exactly solvable, the solution being 
obtained by solving the generalized eigenvalue problem 



QqZv + A J duj'u(uj')^ v (uj) = vZ v 



For v outside the interval lv 2 varies, we get 



Thus for nontrivial solutions 




Nonnegativity of the Hamiltonian excludes negative values for v which lead to 
instabilities, hence we must impose the constraint (cf. [15] eq.(24)) 



n 2 >\ 2 [ dcv^± (2.13) 
J u 

which inter alia implies that at the origin u 2 (u>) is o(u> 2 ) and if the range 
of uj is infinite u 2 (uj) cannot grow at infinity more rapidly than | oj | 2_e for 
< e < 1. Clearly these considirations imply that Q,q cannot vanish, otherwise 
the model will necessarily lead to instabilities. Note however that, had we taken 
Hi = J2 k ut(qk — q) 2 the problem would not arise as should be replaced by 

Although the model is solvable in principle, we will consider weak coupling 
of E with R } and we will obtain the kinetic equation for the oscillator by 
employing the general formalism of [5] that was briefly discussed in section 4 
of paper I and Appendix 3 here, that is, compute the generator $ of Appendix 



3. Because of the separability of Hi, eq.(2.11b), we can use the corresponding 
general formulas of [5] eqs (4.22), (4.22'). To this end we need 

q(t) - e tLst q = q cos n t + ^- sin n t (2.14a) 

W(t) = e tLRt W = J duo u(uj)(t){uj,t) = 

dtu u(oj)((j)(uj) cos ujt-\ smut) (2.146) 

where we write = i{H%,-} etc for the various Liouville operators 

corresponding to (2.1) and we used that for A = 0, (2.12) gives the solution 
q(t),(f)(uj,t). In the present case eqs (4.22), (4.22') of [5] for the probability 
distribution / of the harmonic oscillator are (cf. (2.11b)): 

% = {Hz~ A 2 F, f}+ 
+ T E^-H( P ^)*' { p »«> f}}+~9M{(Pnq)*, (P n q)f}) (2.15) 



2 

n 



F=^^^oj n ){(P n q)*,P n q}+g(oj n )(P n q)*P n q) (2.16) 



n 



Here P n ,iv n are the eigenprojections and eigenvalues of L^, 



/+oo |i 
dse iujs h(s), h(s) =< W W{s) > (2.17a)| 
-oo II 

/+oo 
ds e^ s g(s), g(s) =< W W(s) > (2.176) 
-oo 

r+oo 

Ji{uj)= / dse luJs h(s) (2.17c) 



o 



g(co) = I ds e lujs g(s) (2.17d) 



o 



and < > denotes the average over the chain, assumed in a canonical 
equilibrium state 

p R = I e -4/«M*V)WV(uO)) = L e -PH* (2.18a) 



Z = j <50<5tt e~P HR 



clearly 



u) n = nflo, n 6 Z 

P n A(9,() = ^f d6' e^ e -^ A(6' X) 
with (,0 the action-angle variables of the oscillator: 



(2.186) 



(2.19) 



A simple calculation using (2.19), (2.20) yields 



(2.20) 



(2.21) 



On the other hand explicit calculation of (2.17) involves certain simple 
functional integrals and derivatives. The details are given in Appendix 1. The 
results are: 



~ 2ttu 2 (uj) ~ 

h{v) = 2 , g{u) = i(3u)h{u)) 



(2.22) 



T . , h(uj) . f , ,u 2 (tv') oj 



2 

9{u) 



j3u)' u) 2 — to 



12 



- \ dw' 



to 2 — to' 2 



(2.23a) 



(2.236) 



Notice that because of (2.18), the second of (2.22) follows from the first and 
proposition 4.5 of [5]. An elementary calculation using (2.21), (2.22) gives 

dissipative part of (2.15)= 



A Wqyi r i d 2 f d 2 f d a ' 

= ww i w ^ 

Similarly (2.21), (2.23) substituted in (2.16) give 

u 2 (iv) 1 



(2.24) 



u) 2 n 2 -iv 2 



= constant 



1 -J2gM(P nq yp nq = -H^ f du>^ ] 

n J 



UJ 2 



Since the first does not contribute to the first term on the r.h.s. of (2.15) we 
may write (2.16) as: 

F= MM Hs , A(no) = J- fto-f^ (2.25) 

where we used that the range of u> is symmetric about the origin. By (2.24), 
(2.25), equation (2.15) takes the form 



df 

dt 



\ 2 iru 2 (fl ) ( 8 ( 1 df nn \ d (df nn ,\\ /nnn . 

+ wtUw +w ) + ^ +w )) (2 - 26) 

we may notice that as expected, the Maxwell-Boltzman (MB) distribution 
e -/3#£ j s a stationary solution of (2.26) and F is an integral of the unperturbed 
motion, results that are special cases of proposition 4.1, 4.2 of the general 
formalism of [5] that led to (2.15). 

It may also be remarked that in the litterature, a different kinetic equation 
has been derived for this model ([21] eq.(20) in connection with eq(31) of [15]) 

df a^o) 2 df df 



2 d fn u 2 (n ) r df x(Q ) df 

u 2 (iv) 



U) 2 (u) — Qq) 

In fact it can be shown that this equation is the second order term in the A- 
expansion of the so-called Generalized Master equation (GME), that follows 
by projecting the state of the total system E + R to that of E, namely eqs. 
(3.7') or (5.1) of [5] (see also Appendix 3). Clearly this equation has not 
a nonegative-definite 2nd order coefficient matrix and therefore it does not 
conserve the positivity of / and for that matter obeys no H-theorem. Moreover 
it does not have the MB distribution ce _/3jffs as an equilibrium solution (cf. [5] 
sections 4, 5, [23] section 4, 5). A somewhat different equation but still with 
the structure of (2.27), hence the same undesirable features, has been obtained 
by a similar analysis of the GME, ([30], [31]; see also [34] and the discussion 



in [5] section [5]). Sometimes the mixed derivatives-term is neglected arguing 
that it is "small" ([21] p. 60) so that a conventional Fokker-Planck equation 
results. However even in this case the MB distribution is not stationary, that 
is even for a chain at canonical equilibrium, the external oscillator does not 
evolve toward such an equilibrium state! Moreover it should be noticed that in 
this case this neglect leads to different coefficients for the conventional part of 
the Fokker-Plank equation (compare (2.26), (2.27)). This is a general feature 
that can be understood on the basis of the general formalism presented in [5], 
since the markovian approximation of the GME leads to an equation ((3.7') of 
[5] for the generator O in Appendix 3) totally different from that obtained by 
our formalism ((3.11') of [5] for the generator $ in Appendix 3). This fact has 
been verified in other models as well (e.g. compare eqs (4.1), (4.2) with (4.4), 
(4.5) of [33] section 4). A further argument in favor of (2.26) is provided in 
the next section, where we show that (2.26) is the unique classical limit of the 
kinetic equation of the Lindblad type for the corresponding quantum system, 
which is well-known and widely used, e.g. in quantum optics. 



3. KINETIC EQUATION FOR A QUANTUM HARMONIC 
OSCILLATOR, WEAKLY-COUPLED TO A HARMONIC 
CHAIN 

In this section we consider the quantum system corresponding to the 
classical model of the previous section. It is easily seen that the Hamiltonian 
has the form (hated quantities denoting operators-cf. paper I section 2 for the 
notation) : 

H = Hx + H R + \H I = 
= hQ a + a+'^hLVka^ak + X^2h(eka^ +e* k a k ){a+a + ) (3.1) 

k k 

where 



a = i ( ^ +i i a+ = i ( ^-^' (3 ' 2) 

are the creation-annihilation operators for the oscillator and a k ,cik the 
corresponding quantities for the k-th. bath oscillator, defined similarly^ 1 ), and 
Oq , Wfc are nonegative. 



( 1 ) Non-unit masses can be absorbed into q, p by the same rescaling we did for the classical 
system. 



In fact it is readily seen that any generalized Wigner transformation 
sending functions of g or p to functions of q or p respectively (i.e. (3.17) 
of paper I holds) maps (3.1) to (2.4) (see the comments following it), provided 
6^ is real. For an infinite chain, in the thermodynamic limit we make the 
replacements 



U>k >W Cfc 



— >e(cu) ^ — ► J dw<j(oo) 



ct(uj) being the spectral density of the chain. Consequently, if a = 1, i.e. 
e(k) = e(ujk), then (3.1) is mapped to (2.11) with u 2 (uj) <-» 4e 2 (uj)u}ftQ. A 
detailed comparison of the results of the two sections obviously requires on > 
in the previous section and is discussed at the end of Appendix 1. 

In the litterature, especially in quantum optics, the Hamiltonian (3.1) is 
often treated in the so-called rotating-wave approximation where the terms 
a^a + , ata are neglected (see e.g. [8] p. 336 and eq. (6.2.39), [9] ch. 3 
section 3b, implicity in [12] ch. 12, [22], [36] section 5.14). Sometimes such 
an approximation is made at the end of the derivation of the kinetic equation 
([11] p. 743 eq. (3.1)). To derive a Markovian master equation for weakly 
coupled systems, such an approximation is not necessary and for our purpose 
not desirable, since the Hamiltonian should reduce to the classical one by a 
generalized Wigner transformation (see also Appendix 3). We indicate here 
the steps to apply the formalism of [5] that is, compute the generator $ in 
Appendix 3, noting that H% is separable (cf. [35] for rigorous estimates of the 
approach to equilibrium). Specifically we apply (5.8) of paper I to (3.1), with 

Hj = Wq = J2 WkQ = J2 + ) (fi + « + ) (3-3) 



if v = niVoj^n are the eigenvalues and eigenprojections of j-[Hz,-] = £e , 
n G Z, then clearly 



where P n are the eigenprojections of the number operator N = a + a. Using the 
well-known expression for its eigenvectors 

N\n >= n\n > (3.4a) 

a\n >= \fn \n — 1 > a + \n >= \Jn + 1 |n + 1 > (3.46) 



a simple calculation gives 

T n q = (aSn-! + a + 5 n ,i) (3.5) 
On the other hand, similar expressions to (3.4) for the bath oscillator and 



. Hi 



e iL R s A = ei =^s Ae -i^s } Lr = ^g R , .] 

give 

e iLRS a k = e~ luJkS a k e lLRS a+ = e tuJkS a+ (3.6)| 

Therefore by (3.3) 

W k (s) = e iL * s W k = h(e k e i ^ s a+ + e* k e~ iuikS a+) (3.7) 
As in the previous section, we assume the bath to be in canonical equilibrium 

PR = ^- (3.8) 

A simple calculation using (3.1) gives 

Tr(e-P &R ) = Y[Tr(e-^ k ^ &k ) = ]J (1 - e'^)' 1 
k k 



< a k a k ' >=< a k a k , >= (3.9a) 

_i_ S k k ' n k 

< ala k , >= = — 6 k>k , (3.%) 



< a fc /aj >= \^Y + ljS k>k > (3.9c) 

where < A >= Tr(pRA) and the ^-dependence in (3.9) is shown explicitly so 
that n k has a finite value (PuJk)" 1 in the classical limit. Using (3.9) we readily 
find 

h(s) =< W+W(s) >= h 2 £ |e fc | 2 (e^ kS ^ + l) + e^ fcS y ) 
Using (A. 1.3) we get 

+oo 

~h(co) = J ds e lu}S h{s) = 



27T/1 2 J2 M 2 ( + l) + <"fc) + y - w fc )) (3.10a) 



+00 



u)k Ti u) — u)k 



(3.106) 



We are now ready to obtain the kinetic equation for the density matrix p of the 
oscillator, by applying the above results to the general kinetic equation (5.8) 
of paper I, which for the present model reads: 



dp l r - ^ iX 2 



+ 



(3.11) 



In the thermodynamic limit of an infinite chain, applying (3.10), (3.5) to (3.11) 
using that on > and making straightforward reductions we get ([26] ch. Ill) 



2 A(^o) 



dp 
dt 



fir 



[#E,p] + 



+7TA 2 |e(^ )| 2 a(^o) [^j^- ([a + p, a] + [a+,pa]) + 

(^ + l)([ap,a+) + [a,pa+) 

+00 

A(J2 )= / dujf—^—- + W)| 2 *M 

J \iO — ILq CO + liQ / 



(3.12) 



(3.13) 



As already mentioned, (3.13) is identical with the well-known kinetic 
equation used in quantum optics and derived by making the rotating-wave 
approximation ([8] eq.(6.2.59), [9] eq.(3b.7) p. 118, [11] eq.(3.1), [36] eq. 
(5.1.59)). On the other hand (3.13) is of the form (5.8) of paper I, with the 
following identifications: 



h a p(u>) <-> S n , 



m 1 



n, m = 1, 2 



a % n(Q ) A + 



Bo 



n7 v — + 1 a 



(3.14a) 



7 2 = 7r|e(ft )| 2 ^o) (3-146) 

Therefore its phase-space representation via any involutive generalized Wigner 
transformation (cf.(3.16) of paper I) follows from (5.17) of paper I. From (3.2) 
we find the Wigner transforms of B\,B2 



and consequently their Fourier transforms (paper I, section2) 



Therefore (3.14) above and (5.18) of paper I finally give 

B(a',a) = -(noA*(0*(0*'fa)W) + ^/(^'(OS(v)S(v') 

+iD 2 (sWWiriMrf) - S'iWMSWirfj) (3.15) 

D 1 = 27T 2 7 2 (2n(^ ) + ^) , D 2 = It^^U (3.16) 

Substitution of (3.15) in (5.17) of paper I gives after a straightforward 
calculation and with the aid of (2.12) of paper I, the phase-space representation 
of (3.12): 



fir 



+*D 2 | d(7e <«p(a)(2/i + t/i( W +pe)-^^^-/^^^)) (3-17) 



where [i = ^? [ , ] is the bracket with respect to the *q- product, Hy, is the 
classical Hamiltonian of the oscillator and the kernel of the generalized Wigner 
transformation is 

Q(a) = e x{a) (3.18) 



since ft is an entire (analytic) function without zeros (cf. paper I, section 
2). Eq.(3.18) follows then from the Weierstrass factor theorem for several 
variables (see e.g. [27]). Since the generalized Wigner transformation is 
assumed involutive, so that 

n(a) = n*(-a) & xHr) = x » 

(cf. (3.16) of paper I), we get 



d{—rf) V drj 

that is the Fourier transform of p(a) (and for that matter, of p(a) dX Q^ ) 
is an imaginary function hence it does not contribute to (3.17). The same is 
true for the second term in the first two integrals. Using this, simple reductions 
transform (3.17) with the aid of (3.16), to 

(i lo (n( S l ) + -)- +P p +P -ep) + 



d ( 1 / ,^ , h\dp * \ 



(3.19) 



where 



= ^J e taz x(°)d* (3.20) 

and is the convolution product. This is the desired phase-space equation 
for the quantum oscillator, obtained via any involutive generalized Wigner 
transformation. It is now trivial to obtain the classical limit of (3.19) (or 
equivalently of (3.12)), since all quantities, except \I/ are /i-independent. 
However assuming that 4j , ] is a deformation of the Poisson bracket, 
proposition 3.2 of paper I and (3.18) imply that if (q,p) is mapped to (q,p) 
then lim ^(z) = 0, hence if lim p = pq, then 

7i->0+ R— 0+ 



dt 



(3 - 2i; 



Taking account of (3.9b), (3.14d) and the remarks at the end of Appendix 1 we 
see that (3.21) with a = 1 is identical with the classical equation (2.26) since 
u 2 (u>) } u) > 0, corresponds to 4e(uj)uj£lo as mentioned at the beginning of this 
section. When a ^ 1 the comparison should be made with the kinetic equation 
following by using the Hamiltonian (2.8) instead of (2.11) in section 2. However 
in view of the comments made at the end of Appendix 1, once again the two 
equations turn to be identical. 



4. KINETIC EQUATION FOR MORE GENERAL HARMONIC 
OSCILLATOR MODELS 

In the previous section we obtained the kinetic equation for a quantum 
harmonic oscillator, weakly coupled to an equilibrium bath of other oscillators, 
eq(3.12) within the context of the general formalism of [5]. Moreover its 
phase-space repsesentation via an arbitrary involutive generalized Wigner 
transformation was obtained, eq.(3.19), from which the uniqueness of its 
classical limit, eq. (3.21), follows. Direct application of (3.2) to (3.12) shows 
that the latter is a special case of the following more general equation 



This equation has been considered in the litterature as a master equation 
in quantum optics, describing an electromagnetic field mode interacting with 
an equilibrium bath of bosons. It also includes many other master equations 
used to study open systems in heavy ion collisions ([13] section 3, [14] section 
1,2 and references therein). As it will be seen, it is of the Lindblad type. In 
fact, it follows from it (eq.(5.10")) of paper I) assuming that as functions of 
q,p, its non dissipative part is quadratic and the operators appearing in the 
dissipative part, are linear ([13] section 3, eq.(3.6)). 

In [14] its phase-space representation by means of Wigner's transformation 
or its generalizations corresponding to antinormal and normal ordering of 
operators (Glauber repsesentation) have been given, showing that it is an 
equation of the Fokker-Planck type, i.e. with nonegative-definite 2nd order 
term, provided that 



dp 
dt 



j\H^f>] ~ jt(A + K)[q, [p } p}+] + -(A - K)\p, [p } q}+] 



(4.1) 




£>i,-D 2 >0 D 1 D 2 -D 2 > 



2 A 2 



(4.2) 



4 



The formalism of paper I allows us to generalize this result obtaining the phase- 
space analogue of (4.1) for any involutive generalized Wigner transformation. 
The classical limit of (4.1) is then a simple matter, and as for the special case 
of section 3, it turns out to be unique. 

Rearranging the second term on the r.h.s. of (4.1), we get 

4pi[g, [q,p]]+D 2 \p, \p,p]] - D([q, [p } p]] + \p, [q,p]]) 
n 

-^(foM+l-ftlM+l) (4-1') 

Comparison with (5.5) of paper I shows that (4.1') is of the same form with 
the identifications 

- ( D * ~ D \ ^ -f° A \ 

2 a P~\-D D 2 )' 2 9a(3 ~\-A OJ 

V 1 = V+=q , V 2 = V+=p 

Therefore it can be put in the form (5.8) of paper I, hence it is of the Lindblad 
type, (eq.(5.10") of paper I) with Vi,V 2 as above and 

A ' l 1 ( D * '^-D\ 

V?** = & {-S* - D \ 2 ) < 4 ' 3 > 

provided the latter is nonegative- definite, from which (4-2) follows). Hence its 
phase-space representation is given by (5.17) of paper I with 

B 1 = V 1 = n w 1 (V 1 ) = q, B 2 = V 2 = n~ 1 (V 2 )=p 

and fl^ 1 indicating the Wigner transformation (c.f. section 2.1 of paper I). 
Consequently their Fourier transforms are 

Bifa'.O = 2 7 n5'(7 1 ')5(a B 2 fa,0 = 2^)5' {£) (4.4)| 

so that with the usual notation a = (rj, £), z = (q,p) etc, (4.3), (4.4) substituted 
in (5.18) of paper I, give 

^B(a,a f )= J2 K(a')h nm V m (a) = 



ih A ?h A \ 

+ (— -D)6'{rf)6{r,)6{?)6' {£)-( — + D)6'{r,)6(rf)6{Z)6'(?)) (4.5) 

Substituting (4.5) to (5.17) of paper I gives the desired kinetic equation in 
phase-space. The calculations are tedious but not difficult in principle and are 
given in Appendix 2. The result is 

dp _ jj ® 2 P jj ® 2 P 2D ® 2p 
dt 1 dp 2 2 dq 2 dpdq 

+ |- m 2 q + (A + K )p) (p + ^ p)} + ^ [(- £ + (A - K )q) (p + ^ p)} (4.6) 

where as in the previous section \I> is given by (3.18), (3.20). Therefore it is 
clear that if the generalized Moyal bracket defined by the generalized Wigner 
transformation is a deformation of the Poisson bracket, and (q,p) are mapped 

to (q,p) then, as in the previous section lim\I/(g,p) = and (4.6) has the 

h — >o 

unique classical limit (lim p = po) 

h— >o 

where H% is the classical hamiltonian for the oscillator, (2.4a) 

To recover known results, we explicitly calculate the convolution terms in 
(4.6) for 

Q( v , = e a ^ e+ i ] a e1Z (4.8) 
The following cases are included in (4.8) ([28]) 
a = : Weyl ordering 

a = — 1 : normal ordering (Glauber representation) 
a = 1 : antinormal ordering 
From (4.8), (3.20) we find 

*(«,P) = + n 2 6(q)6"(p)) 



Since x5"{x) = —2S'(x), x5(x) = a direct calculation gives 

| + (A + P)] + + (A - rt] = 

= --«o(A + K )^-^(A- K )^ 
Substituting this in (4.6), we finally get 

-£ = (D 1 - aft K - J n )—^ + (D 2 - ah K oo > )-^ + 2D- 



dt v 2 u/ <9p 2 2Q ' dq 2 dqdp 

+ ^-m 2 q + (A + K)p)p] + |-[(-- + (A - «)g)p] (4.9) 
Introducing the variables 



/n P 



2/T V / 2M^ 

So that x = x\ + 1X2 are the eigenvalues of a corresponding to coherent states 
(e.g. [8] p.107-108), (4.9) becomes 

dp _ ( D 1 {A + k) \ d 2 p f D 2 fl (A - K) \d 2 p D d 2 p 



<9t V 2 ^o 4 J dx 2 \ 2h 4 y<9xf h dx!dx 2 

d d 
-—[((A - k)xi - u;x 2 )p] + — - [((A + k)x 2 + wxi)p] (4.9') 
axi ox 2 

For a = 0, ±1 this is identical with eq. (4.5) and Table III of [14], see also [13] 
eq.(5.19) for the case a = 



In the present paper the general formalism of paper I for the phase-space 
representation of quantum kinetic equations, has been applied to models of 
a harmonic oscillator damped by its interaction with an equilibrium bath of 
other harmonic oscillators. These equations are of the Lindblad type, but 
it should be emphasized that the methods of paper I allow for a phase-space 
representation of any kinetic equation the evolution operator of which is formed 
by algebraic operations between quantum operators. Other applications, as well 
as some more fundamental problems of kinetic theory, already briefly discussed 
in paper I section 1 and [29] section 4, will be examined in the third paper of 
this series. 



APPENDIX 1 



Here we derive eqs (2.22), (2.23) by using simple functional integration 
(see e.g. [24]): By (2.18), Hr is quadratic in 0, tt hence by the formula 



-PH R 



5<j)(uj) V / 5(j)(oS) 
and similarly for ir(uj) we get 

cf)(uj) 5p R 



5<f)5ir (f)(uj)(f)(uj')pR = — J 5(f)57r 



pu' 2 8<t>{u)') 



where we have used integration by parts and have taken into account that 

Similar calculations give 

/5(f)57T tt(uj)7t(uj')pr = -^5{ijJ — u/) 

j 54>Stt ir(uj)(f)(uj')pR = 



1 

Using this and (2.14b) in (2.17a) we find 

his) = ( du) 11 }^} cos a; 5 {A. 1.1) 

J (3u) z 

Similarly, using 

8<t>{u)') 5it(uj') { h 5it(uj') 
and that / ScpSn pr = 1 , we find from (2.17b) that 

g(s)= ! duj^-^- sinus (A. 1.2) 



Using the well-known formulas 

+00 +00 

J ds e ias = 2n6{a), J ds e ias = nS{a) + % - 



we get 



i-00 

J ds e ias cos u)s = ir (s(a - w) + S(a + co)j (Al.3a) 

-00 

+00 

y e ias sina;s = n(s{a-uj) -5{a + u;)^ (A.l.Sb) 

-00 
+00 

/7T / \ ZCl 
ds e ms cos a;s = - [5(a-co)+S(a+co)J-\ — (A. 1.3c) 





+00 



[ dse ias smu;s=—(s(a-uj)-S(a+u;)) (A.1.3d) 

J 2 V / a 2 — lv z 



Substituting (A. 1.1) in (2.17a) and using (A. 1.3a) we get 

h( a ) =tt J ^(a-iv) + 5(a + iv)^^j^div {AAA) 

which gives the first of (2.22). The expressions for g,h,g in (2.22), (2.23') 
are similarly obtained. We may notice here that if instead of the Hamiltonian 
(2.11) we consider (2.8), then nothing changes in the preceding calculations as 
well as those of section 2, except that the integration in (A. 1.1), (A. 1.2) are 
with respect to 6. Hence in (A. 1.3) on is replaced by oo(9) and consequently by 
a change of variables 6 — >u>(0), (A. 1.4) gives 



provided u>{6) is an invertible, even function of 9 (notice that because of (2.10), 
u){6) is necessarily either even or odd). It is now clear that in the kinetic 
equation (2.26) (and in (2.25) as well), u 2 is replaced by u 2 a, which has already 
been remarked in section 2. 



If the spectral variable uj in (2.11) is nonegative then u(uo) can be extended 
by putting u(ui) = 0, for uj < and all calculations remain unaltered except 
that 

(i) in (A. 1.4) (and the corresponding equation for g), the one 5-function does 
not contribute 

(ii) by (2.22), h(uj-i) = g(oo-i) = 0, hence the 2nd term of (2.21) does not 
contribute to (2.15). 

Then it is readily checked that these changes imply a factor | on the r.h.s. of 
(2.24). The case uj > is used in section 3 to show explicitly that (2.26) thus 
modified, is the classical limit of the kinetic equation for the corresponding 
quantum system. 



APPENDIX 2 

Here we derive (4.6). For the dissipative part we substitute (4.5) to 
eq(5.17) of paper I. Specifically we compute 



2A 2 



3*2 



(27T)°ft 



3? / dada'da" 



i(a+a'+a") 

Q(a + a' + a") 



B(a',a) 



p w {(j")e^ a ' ACT ") smhp(a A (a' + a")) 



where p w is the Fourier transform of the Wigner transform of p, a = (//,£), 
a A a' = r/^ — r/£ (see section 2 of paper I for the notation). Substitution of 
B(a',a) from (4.5) and straightforward integrations of the 5-functions give 



271-^ 



ft / dae iaz p(a) 



D\ ( H 2 £ 2 + ip£q - ti£dr,x(<7 ) ) + 



+D 2 {p 2 r] 2 - iprjp + fir]d v x(<rj) ~ ~ d )[p + »WI + P 2f lC ~ PVd v X^)) 



where we used (2.12) of paper I and (3.18). It is not difficult to see that this is 



_ d 2 p _ d' 2 p _ d' 2 p . d r 



dqdp dp 



+A— 

dq 



(A.2.1] 



For the nondissipative term of (4.1') we proceed as follows: Since Hy, is given by 

2 

(3.1), (3.2) its Wigner transform is the classical Hamiltonian = \ + ti^q 2 
hence its Fourier transform is 



6s = -^(IWLM 



Substituting this in (5.12b) of paper I we get after some simple reductions and 
using (2.12) of paper I: 



= h\\ J daeiaZ P( a ){ Ai WP-^r]d< i x{°) / 

+ ^ J dae™ z p{a)^q-^d n x{v)\ 
This is easily transformed to 

-*[*„,,] = { H E , P} + tf * ( * ,) * ( * @P ) (A.2.2) 



Finally, since (g,p) is mapped to (q,p) by the Wigner transformation, a similar 
calculation using (5.15b) of paper I yields 



i 



, , 1 (dp dp d I * ,\ d / * ,\\ , . „ „, 



Combining (A. 2. 1-3) gives (4.6). 



APPENDIX 3 

In sections 2, 3 we derive kinetic equations for a harmonic oscillator, 
weakly-coupled to a harmonic chain, by applying the general formal theory of 
open systems S weakly coupled to equilibrium baths R, developped in [5] (see 
also paper I, section 4.2). The starting point of this formalism is the solution 
of the Liouville (or von Neumann) equation for the total system, for positive 
times. This implies that the corresponding resolvent operator is analytic on 
the complex upper half-plane and has a cut along (part of) the real axis, 
which contains its necessarily existing continuous spectrum. Assuming that 
the restriction of this resolvent to the state space of the open system admits 
a meromorphic analytic continuation through the cut in the lower half-plane, 



it can be shown that the corresponding exponentially decaying contributions 
to the solution of the exact dynamics satisfy a Markovian equation. Moreover, 
its generator can be given a representation in the subspace of the open system 
under certain assumptions which will not be discussed here (cf. [5] section 
2). This generator, $ say, has been computed to second order in the coupling 
parameter of the open system with the bath ([5] section 3) and it has been 
shown to be identical to that proposed by Davies ([3] eq.(2.17). See also [32]). 
It is used in this paper for the particular system studied. Specifically if the 
Liouville operator L of the total system is splitted as (see (4.6) of paper I) 

L = L s + L R + ALj = L + AL 7 

and P is the projection onto the subspace of the open system (see (4.7) of paper 
I), then to second order in the coupling parameter A, $ is ([5] eq.(3.11')) 

T +oo 

$ = -iPLP - A 2 t lim ^ J dt j ds Pe iLo{s+t) L I Qe- iL ° s QL I e' iL ° t P 

-T 

where Q is the complementary projection to P. Here it has been assumed that 
P commutes with Lq, i.e. for A = 0, E evolves independently of R which stays 
in equilibrium, and that Ls has a point spectrum. 

On the other hand, projecting Liouville's (or von Neumann's) equation 
to the P subspace, the so-called Generalized Master Equation (GME) results, 
which under the same assumptions, lead to a quite different generator, say, 
([33] section 2, eq.(2.10) and [5] eq.(3.7')) 



+oo 

9 = -iPLP - A 2 J ds PL I Qe- lLoS QL I e tLoS P 
o 



Notice that 



T 

$= lim — [ dt e iLot Qe- lLat 
T^+oo 2T J 

-T 

is essentially the generator that lead to kinetic equations with serious defects 
(nonconservation of positivity of the states and for that matter, no H-theorem), 
like (2.27) for classical systems (see the discussion in [5] section 5) and which call 
for additional assumptions, like the rotating- wave approximation mentioned in 
section 3, in order to get rid of these defects. 

For a separable interaction Hamiltonian, $ takes the form (4.22') of [5], 
which for the model of section 2 gives (2.15), whereas for its quantum analogue 
it gives (3.11) (notice that PLjP = in our case). 



REFERENCES 



1) C. Tzanakis and A. P. Grecos, "Generalized Moyal structures in phase space, 
kinetic equations and their classical limit. I: General Formalism" (Preprint 
1996) 

2) G. Lindblad, Comm. Math. Physics 48 (1976) 119 

3) E. B. Davies, Comm. Math. Physics 39 (1974) 91 

4) H. Spohn and J. L. Lebowitz, Adv. Chem. Physics 38 (1978) 109 

5) A. P. Grecos and C. Tzanakis, Physica A 151 (1988) 61 

6) C. Tzanakis, Ph. D. Thesis (Universite Libre de Bruxelles 1987) 

7) R. Diimcke and H. Spohn Zeit. Physik B 34 (1979) 419 

8) W. H. Louisell, "Quantum statistical properties of radiation" (Wiley- Inter- 
science, New York 1973) 

9) F. Haake in "Statistical Treatment of open system by Generalized master 
equations" Springer Tracks in Modern Physics Vol 66 (Springer, Berlin 1973) 

10) H. Haken "Light waves, photons, atoms" (North-Holland, Amsterdam 
1981) 

11) G. S. Agarwal, Phys. Rev A 4 (1971) 739 

12) G. S. Agarwal, "Quantum statistical theories of spontaneous emission and 
their relation to other approaches", Springer Tracks in Modern Physics 
Vol 70 (Springer, Berlin 1974) 

13) A. Sandulescu and H. Scutaru, Ann. Physics 173 (1987) 277 

14) A. Isar, W. Scheid and A. Sandulescu, J. Math. Physics 32 (1991) 2128 

15) P. Ullersma, Physica 32 (1966) 27 

16) P. Ullersma, Physica 32 (1966) 74 

17) R. I. Cukier and P. Mazur, Physica 53 (1971) 157 eq.(l) 

18) G. Ford, M. Kac and P. Mazur, J. Math. Phys. 6 (1965) 504, eq.(2) 

19) D. Vitali and P. Grigolini, Phys. Rev. A 39 (1989) 1486 eq.(2.1) 

20) K. Kobayashi and E. I. Takizawa, Mem. Faculty Engineering, Nagoya 
University 20 (1968) 144, eq.(2.1) 

21) P. Ullersma, Physica 32 (1966) 56 

22) H. Carmichael, "An open system approach to quantum optics" 
Springer Lecture Notes in Physics m.18 (1993), p. 13, eq (1-47) 

23) C. Tzanakis, "Linear kinetic equations for finite classical systems: 
Explicit form and its mathematical and physical foundations", Preprint, 
University of Crete (1995) 



24) M. H. Zaidi, Forschr. Phys. 31 (1983) 403 

M. J. Beran, "Statistical continuum theories" (Wiley-interscience, New 
York 1968) 

25) G. Rosen, "Formulation of classical and quantum dynamic theory" , 
(Academic press 1969) 

26) P. Hatjimanolaki, Ph.D. Thesis (Universite Libre de Bruxelles 1988) 

27) K. Knopp, "Theory of functions: Part II", (Dover, New York 1947), ch. I 

28) G. V. Dunne, J. Phys. A 21 (1988) 2321 

29) C. Tzanakis, Physica A 179 (1991) 531 

30) S. A. Adelman, J. Chem. Physics 64 (1976) 124, eq (3.19) 

31) N. G. van Kampen and I. Oppenheim, Physica 138A (1986) 231 eq.(3.5) 

32) E. B. Davies, Math. Annalen, 219 (1976) 147 

E. B. Davies, Ann. Inst. H. Poincare A28 (1978) 91 

33) C. Tzanakis and A. P. Grecos, Physica A 149 (1988) 232 

34) N. Corngold, Phys. Rev. A15 (1977) 2454 
N. Corngold, Phys. Rev. A24 (1981) 656 

35) E. B. Davies, Comm. Math. Physics 33 (1973) 171 

36) C. W. Gardiner, "Quantum Noise", (Springer Berlin 1991) 



